Two-particle spatial correlations in superfluid nuclei 
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We discuss the effect of pairing on two-neutron space correlations in deformed nuclei. The spatial 
correlations are described by the pairing tensor in coordinate space calculated in the HFB approach. 
The calculations are done using the D1S Gogny force. We show that the pairing tensor has a rather 
small extension in the relative coordinate, a feature observed earlier in spherical nuclei. It is pointed 
out that in deformed nuclei the coherence length corresponding to the pairing tensor has a pattern 
similar to what we have found previously in spherical nuclei, i.e., it is maximal in the interior of the 
nucleus and then it is decreasing rather fast in the surface region where it reaches a minimal value 
of about 2 fm. This minimal value of the coherence length in the surface is essentially determined 
by the finite size properties of single-particle states in the vicinity of the chemical potential and 
has little to do with enhanced pairing correlations in the nuclear surface. It is shown that in nuclei 
the coherence length is not a good indicator of the intensity of pairing correlations. This feature is 
contrasted with the situation in infinite matter. 



I. INTRODUCTION 

According to pairing models, in open shell nuclei the 
nucleons with energies close to the Fermi level form cor- 
related Cooper pairs. One of the most obvious mani- 
festation of correlated pairs in nuclei is the large cross 
section for two particle transfer. In the HFB approach, 
commonly employed to treat pairing in nuclei, the pair 
transfer amplitude is approximated by the pairing tensor. 
In coordinate space the pairing tensor for like nucleons is 
defined by 

K(ri*i,r2* 2 ) = (HFB\^{r lSl )^{r 2S2 ) \HFB) (1) 

where \HFB) is the HFB ground state wave function 
while 'f? (fs) is the nucleon field operator. By definition, 
the pairing tensor k (f{s\,f2S2) is the probability ampli- 
tude to find in the ground state of the system two corre- 
lated nucleons with the positions f\ and fi and with the 
spins s\ and si- This is the non-trivial part of the two- 
body correlations which is not contained in the Hartree- 
Fock approximation. 

In spite of many HFB calculations done for about half 
a century, there are only few studies dedicated to the 
non-local spatial properties of pairing tensor in atomic 
nuclei PHI]- One of the most interesting properties of 
the pairing tensor revealed recently is its small exten- 
sion in the relative coordinate r = f\ — rj. Thus in Rcf. 
Q it is shown that the averaged relative distance, com- 
monly called the coherence length, has an unexpected 
small value in the surface of spherical nuclei, of about 



2-3 fm. This value is about two times smaller than the 
lowest coherence length in infinite matter. Similar small 
values of coherence length have been obtained later for 
some spherical nuclei |J| and for a slab of non-uniform 
neutron matter [f|. 

The scope of this paper is to extend the study done in Rcf. 
[|| and to investigate axially-deformed nuclei. It will be 
shown that in axially-deformed nuclei the pairing tensor 
has similar spatial features as in spherical nuclei, includ- 
ing a small coherence length in the nuclear surface. The 
paper is organized as follows. In section [Hi the general 
expression of the pairing tensor is derived in an axially 
deformed harmonic oscillator basis. Expressions of the 
pairing tensor coupled to a total spin S=0 or 1 and as- 
sociated projection are also presented in three particular 
geometrical configurations. In section lHll local as well as 
non-local part of the pairing tensor are discussed for few 
axially deformed nuclei, namely 152 Sm, 102 Sr and 238 U. 
Results concerning coherence length are also presented 
and interpreted in a less exclusive way compared to Ref. 
[H . Summary and conclusions are given in section IIVI 



II. PAIRING TENSOR FOR 
AXIALLY-DEFORMED NUCLEI 

As in Ref.0, we calculate the pairing tensor in the 
HFB approach using the D1S Gogny force @. To de- 
scribe axially-deformed nuclei we take a single-particle 
basis formed by axially-deformed harmonic oscillator 
(HO) wave functions. In this basis the nucleon field op- 
erators can be written as 
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where the HO wave function is 

<f>mu (r) = e me B w „ (?) 



(2) 



(3) 



The quantum numbers m and s are the projection of the 
orbital and spin momenta on symmetry (z) axis; v are 
the radial quantum numbers v = (n±, n z ). The function 

yt\m\v (r) = $l\m\v ( r ±> z ) is gi ven by 

{r±,a x ) (4) 



where 

<Pn z (z,a z ) 
and 
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In the above equations, a z and a± are the HO param- 
eters in the z and perpendicular directions, which are 
related to the HO frequencies by a z = Moj z /K and 
a±_ = Muj^/h, respectively with M the nucleon mass, 
and H nz and L n± are Hermite and Laguerre polynomi- 
als, respectively. 

Using the expansion (2) it can be shown that the pair- 
ing tensor in coordinate representation can be written in 
the following form (the spin up is donoted by "+" and 
the spin down by "-") 



K + , r"a-) = 2J ^|mi|i/i (n) K|mi|iA, (^) 
mi>0 v\ U2 
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In the above expressions we have introduced the pair- 
ing tensor in the HO basis 



ta 1 a 2 — K mivi,m2V2 ~ 2S2(0|c miSl y 1 C- m2 - S2 v 2 |0) 



(9) 



where = mi + si = m 2 + s 2 . 

In the present study we calculate the pairing tensor 
corresponding to three geometrical configurations shown 
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FIG. 1: The geometrical configuration (a) corresponding to 
two neutrons in the xz plane. R and r indicate the c.o.m 
position and the relative distance of the two neutrons. 




FIG. 2: The geometrical configuration (b) corresponding to 
two neutrons in the xy plane. R and r indicate the c.o.m 
position and the relative distance of the two neutrons. 



in Figs{T]|3j they have the advantage of a simple separa- 
tion between the center of mass (c.o.m) R = (rj + r5)/2 
and the relative r '= f\ — f*i coordinates. 

For a finite range force, as the D1S Gogny force used 
here, the pairing tensor has non-zero values for the to- 
tal spin 5 = and 5 = 1. How these two channels are 
related to the pairing tensors <jTj)-<(H]) depends on the ge- 
ometrical configuration. Thus it can be shown that for 
the configuration displayed in FigfTJ the following rela- 
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FIG. 3: The geometrical configuration (c) corresponding to 
two neutrons in the yz plane. R and r indicate the c.o.m 
position and the relative distance of the two neutrons. 
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tions arc satisfied: 

[K(rlsi,rls 2 )] o = V% K(fi+,f2—) (10) 



[K(fis 1 ,f*2S 2 )] 10 = 



(11) 



[K(ns 1 ,i%s 2 )] 11 = [K(fls 1 ,f 2 s 2 )} 1 _ 1 = n(fi+,f 2 '+) 

(12) 

where the notation means that the pairing tensor is 
coupled to total spin S = i with the projection S z = j. 

For the configuration shown in Figj2] the pairing ten- 
sor « (risi,r~2S2) is a complex quantity and we have the 
relations 



K{fiSi,f^s 2 )] 00 = v2 Re(K(fl+,f^-)) 



[k {fisi,fis 2 )} 1Q = iV2 Im(K(fi+,f 2 -)) 



(13) 



(14) 



[K(f{si,f 2 S 2 )} 11 = [n(fisi,f 2 \s 2 )} 1 _ 1 = K{f\ + ,f* 2 +) 

(15) 

Finally, for the configuration (c) of Fig(31 we have 

[K(r!si,rls 2 )] 00 = V2 Re (k (f{+, f 2 -)) (16) 



[ft(rlsi,r1s 2 )] 10 = (17) 



[K(fisi,f 2 \s 2 )} 11 = K(fl+,r" 2 +) (18) 



[n(flsi,r" 2 s 2 )} 1 _ 1 = K(f{-,f* 2 -) = k* (rl+,r1+) (19) 

The results for the pairing tensor shown in this paper 
are obtained by solving the HFB equations in a HO ba- 
sis with 13 major shells for deformed nuclei. We have 
checked that by increasing the dimension of the basis 
the spatial properties of the pairing tensor do not change 
significantly up to distances of about 10 fm in the nu- 
clei studied here. This shows that a finite discrete 13 
major shell HO basis correctly describes these nuclei in 
the domain of interest of this work, in particular that 
continuum coupling effects can be ignored. 



III. RESULTS AND DISCUSSION 
A. Local and non-local parts of the pairing tensor 

We shall start by shortly discussing the local part of 
the neutron pairing tensor. To illustrate the effect of the 
deformation, in FigU]is displayed the neutron local part 
of the pairing tensor for 152 Sm calculated in the spheri- 
cal configuration /? = and in the deformed ground state 
f3 = 0.312, where (3 is the usual dimensionlcss deforma- 
tion parameter. The color scaling on the right side of 
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FIG. 4: (Color online) The local part of the pairing tensor for 
152 Sm. The upper and lower panels correspond to the spheri- 
cal configuration and the deformed ground state, respectively. 



plots indicates the intensity of the local part of the pair- 
ing tensor. In the spherical state the spatial structure of 
the local part of the pairing tensor can be simply traced 
back to the spatial localisation of a few orbitals with 
energies close to the chemical potential 0j- For 152 Sm 
the most important orbitals are 2/ 7 / 2 , 1^9/2, 3p 3 / 2 and 



2/5/ 2 - As seen in FigJH in the deformed state the spatial 
pattern of the pairing tensor is more complicated. This 
stems from the fact that it requires quite many single- 
particle configurations to explain its detailed structure. 
The spatial distribution of the configurations contribut- 
ing the most to the pairing tensor are shown in Fig(S] ; 
the plots correspond to the contribution of single-particle 
states of given f2 and parity, with a different scaling for 
each panel. 

We shall focus now on the spatial structure of the non- 
local neutron pairing tensor. In Figs JB][5] are shown some 
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FIG. 5: (Color online) The spatial structure of the single- 
particle blocks O w which have the largest contribution to the 
pairing tensor shown in the bottom panel of Fig|4] 



typical results of \k (ji, rj | 2 in the three geometrical con- 
figurations (a), (b) and (c) described in FigQ][3]and for 
152 Sm, 102 Sr and 238 U. At the spherical deformation, the 
three geometrical configurations are equivalent. For 102 Sr 

that manifests coexistence feature, |« (Ri^J 1 2 i s shown 

only for the prolate minimum. For 238 U, the ground state 
as well as the isomeric state are displayed. The color scal- 
ing on the right side of plots indicates the intensity of the 
pairing tensor squared multiplied by a factor 10 4 . First 
of all we notice that the pairing tensor for S = 1 (FigJSl 
bottom panel) is much weaker than for S — (FigEl top 
panel) by a factor ~ 20. This is a general feature in open 
shell nuclei (the pairing channel S=l is significant in halo 
nuclei such as 11 Li). Therefore, in what follows we shall 
discuss only the channel 5 = 0. 

Figs 1MH1 show that with deformation the pairing tensor 
is essentially confined along the direction of the c.o.mall 
coordinate. As in spherical nuclei, the pairing tensor can 
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FIG. 6: (Color online) Non local k(R, r) 2 for the isotope 
152 Sm. The deformation is indicated by J3 and S is the spin. 
For the spherical case (upper panel), k(R, r) is averaged over 
the angles of R and r. For the deformed case (lower panel), 
geometrical configuration (a) of FigfT] has been adopted. 



be preferentially concentrated either in the surface or in 
the bulk, depending on the underlying shell structure. 
The most interesting fact seen in Figs. [6][8] is the small 
spreading of pairing tensor in the relative coordinate. 
This is a feature we have already observed in spherical 
nuclei. In Ref. [|| , it is discussed that the predilection for 
small spreading in the relative coordinate is caused by 
parity mixing. We will go back to this point in section 

mm 

Quantitatively the spreading of the pairing tensor in 
the relative coordinate can be measured by the local co- 
herence length (CL) defined in Ref.Q. In the present 
study of deformed nuclei, as particular angular depen- 
dences are assumed according to the three geometrical 
configurations (a) , (b) and (c) , the following formula has 
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FIG. 7: (Color online) Non local K,(R,f) 2 for the isotope 
152 Sin. The deformation is indicated by /9 and S is the spin, 
(a), (b) and (c) indicate the geometrical configurations shown 
in Figs HE] 



FIG. 8: (Color online) Non local n(R,r) 2 for isotopes 102 Sr 
and 238 U. For the latter are shown two cases corresponding 
to the ground state (middle panel) and to the fission isomer 
(bottom panel). Calculations have been made assuming con- 
figuration (a) of Fig[T] 
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been used: 



12r 



'/r 4 |K(i?,r)| 2 dr 
/ t 2 \k{R : r)\ 2 dr 



(20) 



The pairing tensor K(R,r) corresponds to a given total 
spin S= and a given geometrical configuration. For 
spherical nuclear configurations, the expression adopted 
for the CL is the standard one defined in Ref.01 where 
averages are taken over both the angles of R and f. 

In Fig 121 we present the neutron CL calculated for var- 
ious deformed nuclei and configurations (a), (b) and (c) 
described in FigQ][3J We notice that inside the nucleus 
the CL has large values, up to about 10-14 fm. This or- 
der of magnitude was already found in spherical nuclei. 
However, the CL displays much stronger oscillations com- 
pared to spherical nuclei, especially for the geometrical 
configurations (a) and (b). This behaviour can be at- 
tributed to the large number of different orbitals implied 
in pairing properties for deformed nuclei. An interesting 
feature seen in Fig(9]is the pronounced minimum of about 
2 fm far out in the surface which appears for all isotopes 
and all geometrical configurations. The minimum found 
here has a similar magnitude as in spherical nuclei. A 
small coherence length of ~ 2 fm in the surface of nuclei 
we have also found for the protons. In the proton case, 
the Coulomb force has not been taken into account in the 
pairing interaction but it is not expected to change the 
CL strongly. 



B. Discussion of the coherence length 

Compared to the smallest values of the CL in nuclear 
matter, of about 4-5 fm (see Fig[10] for symmetric and 
neutronic matter), the minimal values (~ 2 fm) of the 
CL in nuclei are astonishingly small. The question which 
then arises is what causes such small values of the CL in 
the surface of nuclei. Since, as we have just mentioned, 
the general behaviour of the CL is similar in spherical 
and deformed nuclei, in what follows we shall focus the 
discussion on spherical nuclei. As a benchmark case, we 
will consider the isotope 120 Sn. In this case, the CL will 
be calculated as in Ref.yf, 



Jr 2 \K(R,f)\ 2 d 3 r 
\\ J\K(R,r)\ 2 d s r 



(21) 



where averages are taken over both the angles of R and 
r. 

A possibility explaining the small CL in the surface of 
finite nuclei could be, as, e.g., suggested in Q, that pair- 
ing correlations are particularly strong there. Indeed, 
in the surface the neutron Cooper pairs have approxi- 
matively the same size as the deuteron, a bound pair. 
This is a situation similar to the strong coupling regime 
of pairing correlations. However, it is generally believed 
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FIG. 9: Coherence length for various isotopes. f3 denotes the 
deformation while (a),(b),(c) are the geometrical configura- 
tions shown in Figs. HE] 
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FIG. 10: Coherence length in symmetric and neutronic mat- 
ter according to the density normalized to the saturation den- 
sity and calculated with the D1S Gogny force. 

that nuclei are with respect to pairing in the weak cou- 
pling limit 0. In what follows, we shall examine whether 
there exists a correspondence between the magnitude of 
the CL and an enhancement of pairing correlations in the 
nuclear surface. Even though a local view can only give 
an incomplete picture because of fluctuations, a quantity 
which can be used to explore the spatial distribution of 
pairing correlations is the local pairing energy 

E C {R) = - J d 3 rA(R, r)n(R, r), (22) 

where A(i?,r) is the nonlocal pairing field. In practice 
in Eq. (|22[) . we use the angle-averaged quantities. 

The localisation properties of E C (R) can be seen in 
FigJTT] (black line) where we show the results for several 
spherical nuclei. We notice that in the surface region 
where the minimum of the CL is located there is a lo- 
cal maximum of \E C (R)\ present for all nuclei considered. 
The largest value of E C (R) (in absolute value) is not nec- 
essarily located in the surface region and the oscillations 
of the inner part of the distributions seem mostly due to 
shell fluctuations. 

In order to better exhibit a surface enhancement of 
pairing correlations, we have to consider a normalized 
pairing energy, otherwise the strong fall off of the density 
will mask to a great deal the local increase of pairing. 
One could divide E C (R) by the local density, as done 
in [l[. However, here we prefer to divide by the local 
pairing density k(R) ~ k(R,0), leading to the following 
definition of an average local pairing field 

A av (R) = -J— f d 3 rA(R,r)n(R,r) (23) 
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FIG. 11: (Color online) Pairing correlation energies (right 
scale in MeV) and the average Dairine fields fleft scale in 
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In practice in Eg. ([23)) wc again use the angle-averaged 
quantities. We remark that with a zero range pairing 
force the above definition of the average pairing field 
gives the local pairing field. The localisation properties 
of A av (R) can be seen in FigQl] (red line). We notice 
a qualitatively similar behavior as for E C (R). However, 
due to the normalization, the average pairing field has 
significant values out in the surface region. A closer 
inspection of FigfTT] shows that the averaged pairing 
field reaches by 20% further out into the surface of the 
nucleus compared to the particle density (blue line). 
This can be quantified by the corresponding root mean 
squared values. This push of pairing correlations to 
the external region is determined by the localization 
properties of orbitals from the valence shell, which 
give the main contribution to the pairing tensor and 
pairing field. Since these states are less bound they 
are more spatially extended than the majority of states 
which determine the particle density and the nuclear 
radius. Moreover, the increase of the effective mass in 
the surface also probably plays an important role. Like 
the local pairing energy E C (R), the average pairing field 
A a „(i?) presents a generic local maximum in the surface 
region with a local enhancement of pairing correlations 
(at tenth the matter density in 120 Sn, the average pairing 
field still reaches a relatively large value of ~ 0.5MeV). 
On the other side, this maximum is not necessarily an 
absolute one and higher pairing field values can appear 
in the interior of nuclei, depending on the underlying 
shell structure. 

In order to understand if this local enhancement of 
pairing correlations is able to explain the minimum value 
of ~ 2 fm of the CL in the surface of finite nuclei, we 
have calculated the CL under the same conditions as be- 
fore but with a variable factor a in front of the (S=0, 
T=l) pairing intensity of the D1S Gogny force (and only 
there). The result is shown in Fig|12]for a between 1.0 
and 0.5 (top panel) for 120 Sn. It should be mentioned 
that for a = 1.0, the 120 Sn pairing energy is equal to 
~ 19MeV whereas for a = 0.5 it is ~ 0.5MeV which can 
be considered as a very weak pairing regime. In spite of 
these extreme variations of the pairing field, the values of 
the CL are changing overall very little, except for R < 1 
fm. At R ~ 6 fm, the variation is less than 0.2 fm. As 
we will see, this behavior is completely different in nu- 
clear matter. From this study, it becomes clear that the 
CL is practically independent of the pairing intensity, in 
particular in the surface of finite nuclei. 

Therefore, we must revisit the interpretation proposed 
in our preceding paper 0], that the minimal size of ~ 2 
fm of Cooper pairs in the nuclear surface is a consequence 
of particularly strong local pairing correlations. From the 
fact that a completely different behavior is obtained in 
infinite nuclear matter (see below and Fig fl6|) . the small 
size of the CL in the surface of nuclei seems to be strongly 
related to the finite size of the nucleus. 

At this stage of our analysis, it is important to clar- 
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FIG. 12: (Color online) Coherence length calculated with 
total pairing tensor k (top panel), even n e and odd k part of 
the total pairing tensor (bottom panel) for different intensity 
of pairing strength, in the case of 120 Sn. 



ify the role of parity mixing which was put forward in 
our preceding work [|| on the behaviour of the CL. In 
the bottom panel of Figfl2l is displayed the CL calcu- 
lated either with the even part of the pairing tensor K e 
or with the odd one k , for the same values of a as be- 
fore. One sees that, in both cases, the value of the CL 
does not depend much on the intensity of the pairing. 
This conclusion holds here for all the values of R. Com- 
paring the curves in the two panels of Fig |12l one sees 
that the even/odd CL's have sensibly larger values in the 
center of the nucleus (around ~ 10 fm) than the CL cal- 
culated with the full k (6-8 fm), almost independently of 
the value of a. In the surface region they are practically 
of the same magnitude (2-3 fm) . These results indicate 
that the parity mixing discussed in Ref. @ influences 
the CL essentially for small values of R. Therefore, par- 
ity mixing cannot be the main reason for the small value 
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FIG. 13: (Color online) Non-local part of the pairing tensor 
k 2 , even k 2 and odd k„ part of the non-local part of the pairing 
tensor and the interference term 2k £ k for 120 Sn. 



of the CL in the surface region. 

The trends observed in FigfT2] can be traced back to 
the variations of n 2 as well CIS /^g f 

and the interference 
term 2n e K plotted in FigfT3] One sees that the interfer- 
ence term is large only along the axes r = and R = 0. 
However, in calculating the CL, \k\ 2 is multiplied by a 
factor r 4 in the numerator of Eq. (|2"Tj) . Hence, the large 
values of this interference term near r = axis will not 
come into play significantly. Therefore, as observed pre- 
viously, parity mixing will be significant essentially for 
R<1 fm. 

This observation is confirmed by looking at the quan- 
tity 



X(R,r) 



r 4 | K (fl,r) p 



(24) 



where N(R) = J °° dr r 2 n(R,r) 2 . This quantity, once 
integrated over r yields the square of the CL, Eq. ([2"T]) 
namely, £ 2 (i?) = f™ X(R,r)dr. 

X(R,r) is presented in Figfl4] for four values of R 
namely, 0, 3, 6 and 9 fm corresponding to the interior 
of the nucleus and the vicinity of the surface. The re- 
sults are displayed for various values of the pairing factor 
a. 

Except for R = 0, X(R, r) and hence £(i?) is not really 
sensitive to the strength of the pairing interaction. The 
large dependence of X(R,r) on the pairing strength at 
R = comes from the comparatively large parity mixing 
already mentioned in connection with Fig |131 which is 
negative and maximum in absolute value for r ~ 10 fm. 
Since the parity mixing tends to disappear as the pairing 



strength decreases, the height of the peak at r = 10 fm 
increases. In contrast, for R = 3,6,9 fm the influence 
of the parity mixing is very modest and the behaviour 
of X(R, r) is determined essentially by k 2 and k 2 . From 
R ~ 3 fm to R ~ 6 fm, one observes a sensitive reduc- 
tion of the magnitude of X(R,r) leading to a lowering 
of the CL. In the vicinity of the surface (R > 6 fm), the 
oscillatory behaviour of X(R, r) disappears. Here, single 
particle wave functions have almost reached their expo- 
nential regime. This explains why at R > 6 fm, X(R,r) 
is characterized by only one major peak. The width of 
this major peak is minimum at the nuclear surface. Its 
broadening for R = 9 fm explains the increase of the CL 
beyond the nuclear surface. 

A more global way to analyze the behavior of the CL 
is to consider directly the dependence on R of the nu- 
merator and the denominator of Eq . ([21"]) . This is shown 
in FigQjj] One sees, that independently of the value of 
a (the color code is the same as for Fig fT2|) . the denomi- 
nator decreases faster than the numerator around R = 6 
fm and beyond. This sudden change in the slope of the 
denominator is accountable for the minimum value of CL. 

A similar analysis of the CL and of the influence of 
pairing correlations has been carried out in infinite mat- 
ter. In Fig 021 we show the CL in infinite symmetric 
nuclear matter as a function of the density p normalized 
to its saturation value po , for the same a values as in the 
HFB calculations for finite nuclei. In the nuclear matter 
case, we see that the CL depends very strongly on the 
pairing intensity, whatever the density. For instance, the 
minimum value of CL increases a lot as pairing decreases. 
This behavior can be understood from an approximate 
analytic evaluation of the CL in infinite nuclear matter 
based on the definition Eq. ([2T|) which differs only slightly 
from the usual Pippard expression (lp| (see Appendix [A"| : 



(25) 

where a = \Af\/\cf\, b = kpAp/Ap with Ap and A F 
the pairing field and its derivative for the Fermi momen- 
tum kp. As discussed in Appendix [XJ the correction 
terms in Eq. (l25l) are very small. 

We see that the CL in infinite matter varies approxi- 
matively inversely proportional to the gap at the Fermi 
surface. This behavior is at variance with the results 
in finite nuclei and particularly where the CL shows the 
minimum, see Fig[12] This clearly indicates that the be- 
havior of the CL, in particular the small value obtained 
in the surface of finite nuclei, is strongly influenced by 
the structure of the orbitals and that pairing plays a sec- 
ondary role. 

In order to examine this question in more detail, we 
show in Fig ll7l the extension of completely uncorrclatcd 
pairs made of Hartree-Fock neutron single particle wave 
functions. We use a definition of the extension of the pair 
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FIG. 14: (Color online) X(R, r) for R=0, 3, 6 and 9 fm in the 
case of 120 Sn. 
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FIG. 15: Evolution of numerator and denominator of (,(R) for 
various values of a, in the case of 120 Sn. 



size similar to the CL of Eq. ([2~Tj) namely, 



£orb(R) = 



fr 2 \Ai(R,r)\ 2 d 3 r 



1/2 



1/2 



(26) 



(/|A ? ;(i?,r)|2d3r) 

The imcorrclatcd pair wave function Ai(R, r) is defined 

as 

MR,r) = i-(2* + 1) CZ Uji C^ 

x J2(-) l {2l+ 9 } )1/2 u nl (r/V2)u Nl (V2R) (2?) 



nNl 



x Pi(cos8)(nlNl;0\n a kn/3liO} 

where (7™* iJ * is the component of the (rnkji) neutron 
single-particle orbital on the HO basis function (n a liji). 
This equation is the same as Eq. (3) of Ref. [3[ with 
the matrix K„^* n „ of the pairing tensor replaced with the 
product of the two C coefficients. 

Since £, rb(R) corresponds to two non-interacting neu- 
trons put into the same orbit and coupled to (L=0, S=0), 
it contains only the correlations induced by the confine- 
ment of the single-particle wave functions. As FigfTTI 
shows, £, or b has a pattern rather similar to the global CL 
displayed in Figs. l9land[T2"l except for the 3si/ 2 orbital. 
Thus, provided this orbital is not strongly populated, a 
change in the relative contributions of the single-particle 
states in the pairing tensor, e.g., induced by varying the 
intensity of pairing correlations, will not cause significant 
modifications in t he g lobal CL. This result has also been 
found by Pastore [ll(. From FigJT7]we see that (except 
for 3si/2), tj,orb(R) exhibits a minimum in the surface of 
the order of ~ 3.5 fm. This is indeed small but still larger 
than the 2.3 fm found with Eq.pijl for a = 1 (or 2.5 fm 
for a = 0.5). The reduction by about 30 percent from 
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FIG. 16: (Color online) Coherence length calculated with dif- 
ferent intensity of pairing strength in symmetric nuclear mat- FIG. 17: Coherence length for Hartree-Fock single particle 
ter. orbitals of the neutron valence shell of 120 Sn. 



3.5 fm to 2.3 fm of the minimum of the CL is very likely 
due to the fact that even for very small pairing some 
orbit mixing takes place (remember that the influence 
of pairing is compensated in the ratio of numerator and 
denominator and that the chemical potential becomes 
not necessarily locked to a definite level but may stay 
in-between the levels) . The cross terms of the wave func- 
tions can be negative yielding a possible explanation of 
the effect. Let us also point out that the CL implying 
only the even part of the pairing tensor (or the odd one) , 
see [H, is of the order of ~ 2.7 fm for 120 Sn, see Fig. 12. 
Therefore, there should exist a slight influence of parity 
mixing in the CL calculated with the full n. 

Nonetheless, the above discussion clearly indicates that 
the small value of the CL in the surface of finite nuclei 
is essentially due to the structure of the single particle 
wave functions. Our conclusion is somewhat different 
from the one put forward in our early paper Q. There, 
we had not explored the behavior of the CL as a function 
of the pairing strength, which led us to conclude that 
the small size of Cooper pairs stems from a local strong 
coupling pairing regime. However, the other results and 
conclusions of Ref. |3| still hold. 

One may speculate about the reason for this radically 
different behavior of the CL in nuclei and infinite nu- 
clear matter. One issue which certainly can be invoked, 
is that in macroscopic systems the number of single par- 
ticle states in an energy range of the order of the gap is 
huge whereas in nuclei we only have a few states/MeV. 
In order to examine more precisely such an effect, let 
us consider, for convenience, the example of a spherical 
harmonic oscillator potential. We want to keep the es- 
sential finite size effects but eliminate unessential shell 
effects. It is well known that this can be achieved via the 
so-called Strutinsky smoothing. Single shells are washed 



out and what remains is a continuum model with energy 
as variable, instead of individual discrete quantum states. 
We therefore can write for the pairing tensor in Wigncr 
space: 



n(R,p) 



dEn(E)f(E; R,p) 



(28) 



where k(E) = ueve is the Strutinsky averaged pair- 
ing tensor [l2j and f(E;R,p) is the Strutinsky averaged 
Wigner transform of the density matrix on the energy 
shell E [l3|. Integration of this quantity over energy up 
to the Fermi level yields the Strutinsky averaged density 
matrix in Wigner space. The latter quantity is shown in 
1 of Ref. 
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Fig. 

A particularity of the Strutinsky smoothed spherical 
harmonic oscillator is that all quantities depend on R 
and p only via the classical Hamiltonian H c \,(R,p). Wc 
see that the Wigner transform of the density matrix is ap- 
proximately constant for energies below the Fermi energy 
and drops to zero within a width of order fvjj. The corre- 
sponding density matrix on the energy shell can then be 
obtained from the quantity shown in Ref. flij by differ- 
entiation with respect to energy. We, therefore, deduce 
that f(E;H c \ ) is peaked around E ~ H c \ with a width 
of order fiw. 

The above integral over E in n(R,p) is, therefore, a 
convolution of two functions, one of width ~ A and the 
other one of width ~ hui. As long as the gap is smaller 
than Hu, the R and p behavior of k will be dominated 
by f(E;H c \.), i.e. by the oscillator wave functions. This 
is what happens in finite nuclei. On the contrary, in 
infinite matter or in LDA descriptions, f(E;R,p) is a 
(5-function, S(E — H c \.), and then the R, p behavior of 
k(R,p) is entirely determined by the width of k(E), i.e. 
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by the intensity of pairing. 

This interpretation qualitatively explains the very dif- 
ferent behaviors of the CL with respect to the magnitude 
of pairing in finite nuclei and infinite matter. It also ex- 
plains why the value of the CL in the surface of finite 
nuclei can be much smaller than the one calculated in 
infinite nuclear matter at any density. 

More quantitative investigations along this line are in 
preparation [Til ]. 



IV. CONCLUSIONS 

In this paper we have continued our study of the spa- 
tial properties of pairing correlations in finite nuclei. We 
first generalized our previous work Q to deformed nu- 
clei and found that the spatial behaviour of pairing is 
rather similar to the spherical case. This concerns, for 
instance, the remarkably small value of the coherence 
length (CL) (~ 2 fm) in the nuclear surface. More inside 
the nucleus, sometimes more pronounced differences ap- 
pear. We then concentrated on the reason for this strong 
minimal extension of the CL in the surface of nuclei. It 
was found that this feature is practically independent of 
the intensity of pairing and even seems to survive in the 
limit of very small pairing correlations. A detailed anal- 
ysis of the quantities entering the definition of the CL 
indicates that, in finite nuclei, the latter is mainly deter- 
mined by the single particle wave functions, i.e. by finite 
size effects. This eliminates suggestions that the strong 
observed lowering of the CL in the nuclear surface has 
something to do with especially strong pairing correla- 
tions in the surface Q, or in a surface layer, i.e. with a 
2D effect [l(|. A particular situation seems to prevail in 
the two neutron halo state of 11 Li fl7j . 

We also made the same study in infinite nuclear mat- 
ter. We found that in that case the CL strongly depends 
on the gap and an approximate inverse proportionality 
between the gap and the CL could be established. Con- 
cerning the reason why nuclei and infinite matter behave 
so differently with respect to the CL, we put forward the 
fact that the number of levels in the range of the gap 
value is huge in a macroscopic system whereas there are 
only a handfull of levels in finite nuclei. In such situa- 
tions the numerator and denominator in the definition of 
the coherence length have a similar dependence on pair- 
ing and its influence tends to cancel. From this work, 
it appears that the CL may not be a good indicator of 
the spatial structure of pairing correlations in the case 
of nuclei or of other finite systems with a weak coupling 
situation like certain superconducting ultra small metal- 
lic grains [l8j]. This fact should not make us forget that 
on other quantities pairing in nuclei can have, as well 
known, a strong effect. For example the pairing tensor 
itself, as studied in this work, is very sensitive to parity 
mixing, see Fig |131 where a strong redistribution, i.e. a 
concentration of pairing strength along the c.o.m posi- 
tions of the pairs takes place. Such a feature probably is 



responsible for the str ong enhancement of pair transfer 
into superfluid nuclei [19(. This small extension of the 
pairing tensor in the relative coordinate may not only be 
present in the surface but also in the bulk, depending 
somewhat on the shell structure. However, on average 
a generic but moderate enhancement of pairing correla- 
tions (obtained with the D1S Gogny force) is present in 
the nuclear surface, see FigfTTJ Further elaboration of 
these aspects will be given in a forthcoming paper (l5j . 
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Appendix A: Neutron coherence length in infinite 
nuclear matter 



Introducing the Wigncr transform 

K W (R,k) = [ d 3 r n(R,r)e l%p 



(Al) 



of the HFB neutron pairing tensor k(R, r), the coherence 
length (CL) defined by Eq. (|2"Tj) can be rewritten 



jd*k\^ k K W (R,k)\< 

\l f d 3 k\K W (R,k)\ 2 



(A2) 



In infinite nuclear matter, kw is independent of R, de- 
pends on k only through the length k = \k\, and is 
given by K W (k) = A(k)/2E(k), where A(fc) is the HFB 
neutron pairing field and E(k) = (e(fc) — /j) 2 + A(/c) 2 
the neutron quasiparticle energies with e(fc) the single- 
neutron energies and the neutron chemical potential. 
Substituting these expressions into (|A2jl yields £ nm = 
s/N/D with 

, ^ (e(fc) - (A'(fc)(e(fc) - /,) - A(fc)e'(fc)) 2 
[(e(jfe) - M )2 + A(k) 2 ] 3 



D= / k 2 dk 



A(ky 



(e(fc)- M )2+A(fc)2 

(A3) 

The primed quantities are first derivatives. In order to 
be able to express the integrals analytically, we introduce 
the three following approximations 

1. fx ~ e{kp) = &f where fc^ is the neutron Fermi 
momentum, 

2. e(k) ~ h 2 k 2 /(2m*) where m* is the ^-dependent 
neutron effective mass, 

3. in the usual situation of nuclear physics where the 
gap values are much smaller than the Fermi energy, 
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the functions under the two integrals (|A3[) are suffi- 
ciently peaked around k — k F so that one can take 
AO) ~ A(fc F ) = A F and A'(fc) ~ A'(k F ) = A^. 

Using these assumptions and making the change of vari- 
ables k = xkp, expressions (|A3|) become 



N=a 2 k F 



D= a 2 kl 



oo Jl 



x 2 (x 2 ~\) 2 [b(x 2 - 1) - 2x] 



(•'• 2 



[(a; 2 - l) 2 + a 2 
t- 2 

— Trdx 



2l3 



(A4) 



1) 



with a = \A F /e F \, b = kfA' F /A F . Assuming a ^ 0, the 
integrals on the right hand sides of (|A4[) can be calcu- 
lated analytically using contour integration in the com- 
plex plane and the method of residues (more precisely, 
the integrand for N can be broken into an even function 
for which the integration range, as the one for D, can be 
extended from — oo to +oo and integrated by the method 
of residues, and an odd function which is easily integrated 
after the change of variable y = x 2 ). One gets: 



N=a 2 k f 



2tt aY\ 



1 + Vl + 



-i + VT 



b ( Sir „ < . , a 

Ta T + 3cot (fl) -TT^ 



D= -aki 



l + VTT 



where X and Y are functions of a and b given by 



(A5) 



X = 



Y = 



a 2 b 2 (4,a 2 + 5) - 2(1 + q 2 )(5a 2 + 2) 

64a 2 (l + a 2 ) 2 
a 2 6 2 (21a 4 + 35a 2 + 12) + 4(1 + a 2 )(7a 2 + 4) 



128a 4 (l + a 2 ) 2 



(A6) 



Usually, a is much smaller than one, even at small den- 
sities. Expanding the above expressions around a = 0, 
one obtains 



1 + — (3b' - 126 + 4) + 0K) (A7) 



akpy/2 \ 8 



With a = \A F \/(h 2 kp/2m*), the leading term yields 



1 h 2 k F 
2V2m*\A F \ 



(A8) 



This expression is very close to the Pippard approxima- 
tion of the CL Q3 



£ Pippard 



l h 2 k F 
7r ra*|A F | 



(A9) 



the pre-factor being 1/2V2 ~ 1/2.8 instead of 1/tt. Usual 
values of a and b show that the first correction term in 
(|A7j) is very small. For instance, in symmetric nuclear 
matter at one tenth the normal density, one gets k F ~ 



0.6 far 



.2 and b ~ .3 with the Gogny effective 



force, which yields 3.10 -3 for this term. The next terms 
can be shown to be even smaller. Moreover, numerical 
evaluations of the integrands in (|A3j) for the Gogny force 
show that the three above approximations employed for 
deriving (|A5[) . in particular the third one, are extremely 
well justified for densities ranging from zero to twice the 
normal density in symmetric nuclear matter. 
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